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Abstract. Recent advances in lattice field theory, in computer technology and in chiral perturbation theory 
have enabled lattice QCD to emerge as a powerful quantitative tool in understanding hadron structure. I 
describe recent progress in the computation of the nucleon form factors and moments of parton distribution 
functions, before proceeding to describe lattice studies of the Generalized Parton Distributions (GPDs). 
In particular, I show how lattice studies of GPDs contribute to building a three-dimensional picture of the 
proton. I conclude by describing the prospects for studying the structure of resonances from lattice QCD. 

PACS. 12.38.Gc Lattice QCD calculations - 13.40.Gp Electromagnetic form factors - 13.60.Fz Elastic 
and Compton scattering 



' 1 Introduction 

■ The ability for lattice QCD to make increasingly pre- 
' cise, ab initio calculations of hadron structure has been 

driven by advances in three areas: algorithmic advances, 
enabling both chiral symmetry to be more faithfully repre- 
] sented on the lattice and gauge configurations to be gener- 
' ated more efficiently, computational advances represented 
I by the availability of special-purpose and leadership-class 
. computers and by commodity clusters, and theoretical ad- 
' vances enabling extrapolations to be performed from the 
I quark masses employed in the computation to the physical 

■ up- and down-quark masses. In this talk, I will show the 
impact of these advances on hadron structure, and provide 
a roadmap for future calculations. The layout of the rest 

■ of this talk is as follows. In the next section, I will describe 
' how hadron-structure calculations proceed in lattice QCD. 

I will then review benchmark calculations, such as that of 
the nucleon's electromagnetic form factors. Section |4] will 
describe calculations of moments of Generalized Parton 
Discributions (GPDs) and the impact such calculations 
are having on our understanding of hadron structure. I 
will conclude with prospects for future calculations. 



2 Anatomy of a nucleon-structure calculation 

Hadron structure is expressed through quantities as the 
electromagnetic form factors, describing the distribution 
of charge and currents within a hadron, and the polar- 
ized and unpolarized parton structure functions, describ- 
ing how the longitudinal momentum fraction and spin 
is apportioned amongst the constituents. Recently, Gen- 
eralized Parton Distributions have been introduced[Tl[H 
[3], encompassing both of the concepts above, but allow- 



ing a three-dimensional picture of the nucleon to be con- 
structed. 

GPDs are expressed as the matrix elements of light- 
cone correlation functions Orix): 

Or{x)^l £e'^-g-(Z^)rPe^'^^^-v^'""-^^""^g(^) 

(1) 

where n is a light-cone vector, the parallel-transporter is 
necessary to ensure gauge invariance, and the flavour in- 
dices on the quarks are suppressed. The familiar quark 
polarized and unpolarized parton structure functions are 
then forward matrix elements of this operator, with F — 
In or 7/^75 for the unpolarized and polarized distributions, 
respectively; the GPDs correspond to matrix elements in 
the off-forward direction, with different momenta for the 
incoming and outgoing partons. 

The use of a Euclidean lattice precludes the measure- 
ment of these matrix elements directly. Instead with ap- 
peal to the Operator Product Expansion to expand the op- 
erator 0{x) about the light cone, yielding a set of local op- 
erators which can be measured on a Euclidean lattice, and 
furthermore analytically continued to Minkowski space. In 
particular, the moments, with respect to Bjorken-x, of the 
quark structure functions are obtained in terms of the for- 
ward matrix elements 

^{Mi...p„} ^ ^(^5)^{mj75;^2 _oM„}g (2) 

Specifically, for the unpolarized distribution for a nucleon 
carrying momentum p, we have 

|p) ^ f dxx'^-\{x), (3) 
Jo 

where I have suppressed the Lorentz structure. 
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We have seen in the talk of Christof Gattringer how 
the spectrum is obtained in lattice QCD through the mea- 
surement of so-called two-point functions representing the 
Euclidean correlators of two operators that interpolate be- 
tween a state and the vacuum. Nucleon structure is inves- 
tigated through the measurement of three-point functions, 
illustrated in Figure [1] representing the Euclidean corre- 
lation functions of three operators: 

C3pt(isnk,T;p, q) = 
E e-'P-^-'^-^OlJix, Uk)0(y, r) Jt(0)|0), (4) 

where J is an operator that interpolates between the state 
and the vacuum. Note that the disconnected contribution 
in the right-hand panel of Figure [T] is far more computa- 
tionally demanding that the connected piece; for most of 
the following, I will emphasise the calculation of flavour 
non-singlet, or isovector, quantities, for which the discon- 
nected piece vanishes. 

Following [4] , we insert complete sets of states between 
the operators in eqn. 21 the time-slice sum projects onto 
states of definite momentum, and we obtain: 

n=l,... 

(0|J(0)|n,p)(n,p|O(0)|n,p + q)(n,p + g|jt(0)|0}.(5) 

For sufficiently large r, tgnk ~ C^pt is dominated by the 
lowest-lying state, and we can obtain {n — l,p\0\n ~ l,p) 
after elimination of the vacuum-to-state matrix elements 
from the corresponding two-point functions. 

Computations of hadron structure are being performed 
by several groups, using a variety of fermion and gauge- 
field discretisations. The common theme is calculations at 
increasingly light values of the pion mass, enabling lat- 
tice calculations, together with chiral effective theory, to 
provide insight into physics at the physical quark masses; 
recent reviews are contained in ref. [MS]- During the re- 
mainder of this talk, I will in large part quote results from 
the LHP Collaboration. We employ a hybrid approach, 
using lattices generated using 2 + 1 flavours of Asqtad 
quarks, corresponding to degenerate u/d and strange, by 
the MILC Collaboration^?], but with domain- wall fermions 
(DWF) for the valence quarks, computed on HYP-smeared 
lattices [8]. The generation of lattices using Asqtad fermions 
is very computationally efficient, whilst the DWF valence 
quarks have the desired chiral-symmetry properties, avoid- 
ing the much operator mixing and simplifying the match- 
ing of our results to the continuum. Whilst this approach 
violates unitarity at finite lattice spacing, we expect this 
approach to yield the correct theory in the continuum 
limit [9]. In this talk, I will present results with pion masses 
down to around 360 MeV. 



3 Electromagnetic Form Factors 

These describe the distribution of charge and current within 
the nucleon, and remain the subject of intense experimen- 
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Fig. 2. The points show measurements of the isovector Fi 
form factor at five values of the pion mass; the coloured bands 
are uncertainties from a pole-dominance fit to the data. The 
solid line is a parametrisation of the experimental data[ll]. 

tal and theoretical interest. They are the most straight- 
forward quantity to measure in a lattice calculation, cor- 
responding to the matrix element of the electromagnetic 
current = q-f^iQ. Figure [2] shows the form factor Fi{Q'^) 
determined by the LHP CoUaborationfTIT, together with 
a parametrisation of the isovector experimental data[llj. 
The data show the approach to the experimental parametri- 
sation at decreasing values of the pion mass used in the 
lattice calculation. 

The slope of the form factor at = is related to 
the isovector charge radius (r^)^/^. The left-hand panel in 
Figure [3] shows the charge radius obtained from the cal- 
culation above. A naive linear fit in would clear be 
far below the experimental points, illustrating the need 
to correctly describe the non-analytic behaviour in the 
approach to the physical light-quark masses. The curve 
shows the chiral extrapolation of the charge radius, us- 
ing the finite-range regulator [T2j. The right-hand panel 
shows the isovector mean square charge radius obtained by 
the QCDSF/UKQCD Collaborations, using two fiavours 
of Wilson fermions, with pion masses down to around 
350 MeV; the line is the chiral behaviour in the small- 
scale expansion. [13] 

4 Generalized Parton Distributions 

The electromagnetic form factors described above are par- 
ticular cases of Generalized Parton Distributions. There 
are two GPDs corresponding to the choice _r = 7 • n in 
eqn. [l] and E{x,£,,t), and a further two cor- 

responding to the r = J ■ 7175 in eqn. [U H{x,S,,t) and 
E{x,S,,t). The invariants are: 

t = -A^ ^ -{P-P'f 

f = -n • A/2 (6) 

where is a light-cone vector, and P', P are the four- 
momenta of the incoming and outgoing states. Once again, 
we appeal to the operator product expansion to obtain the 
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Fig. 1. The left- and right-hand panels show the connected and disconnected contributions to the three-point function of eqn.lH 
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Fig. 3. The left-hand panel shows the isovector charge radius using a mixed Asqtad/DWF calculation |10j. together with the 
chiral behavious using a finite-range regulator [12]. The right-hand panel shows the mean-square charge radius using two flavours 
of Wilson fermions'14', together with the chiral behaviour using the small-scale expansion [13). 



matrix elements of local operators, and both the form fac- 
tors described above and the famihar parton distributions 
emerge as special cases: 



and 



ff(x,0,0) = qix) 
H{x,0,0) = Aq{x), 

J dxH{x,£„t) = Fi{t). 



(7) 
(8) 

(9) 



Note here that the x runs from -1 to 1, corresponding to 
antiquark and quark momentum fraction. 

Before proceeding to discuss the new insights enabled 
by the study of GPDs, we begin by discussing the polar- 
ized and unpolarized distributions as a means of bench- 
marking our calculation. A fundamental measure of QCD 
is the nucleon's axial- vector charge gA, corresponding to 
the choice -T = 75; it has additional importance in its role 
as a fundamental low-energy constant of the theory that 
appears in the chiral expansion of other quantities. Fig- 
ure |4] shows a calculation of the axial- vector charge by the 
LHP collaboration, together with the chiral extrapolation 
to the physical quark masses [TS]. The consistency between 
the extrapolated lattice data and the experimental value 



is very encouraging, though the chiral behaviour for this 
quantity is very mild. 

The non-analytic behaviour with decreasing quark mass 
is more strongly exhibited in the unpolarized distribu- 
tions, and in particular in the calculation of the momen- 
tum fraction carried by the valence quarks in the nucleon. 
An analysis of the forward matrix elements for both the 
unpolarized and polarized distributions, as well as for the 
transversity, has been performed using a self-consistent 
chiral expansion [1^; data for the non-singlet momentum 
fraction, and a summary of the results for a range of 
benchmark quantities, is shown in Figure [5l 

4.1 New insights: origin of nucleon spin 

The excitement and signficance of GPDs arises from the 
new insights they can provide into hadron structure, spurring 
new experimental initiatives such as the 12 GeV upgrade 
at Jefferson Laboratory. There is thus a rigorous effort 
both by the LHP Collaboration in a hybrid DWF/Asqtad 
approach [MH, and by the QCDSF/UKQCD Collabo- 
ration using two flavours of dynamical, improved Wilson 
quarks [21] . to measure moments of GPDs. 

An salient example of the new insights facilitated through 
the study of GPDs is how the spin of the nucleon is dis- 
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Fig. 5. The left-hand figure shows the momentum fraction carried by the valence quarks in a nucleon, together with the chiral 
extrapolation to the physical quark masses. The right-hand figure is a summary shows the lattice values at the physical pion 
masses for a range of benchmark quantities, normalized to the experimental value where known (blue) [19]; the errors on the 
corresponding experimental measurements are shown in red. 
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Fig. 4. The red squares show LHPC/MILC data for the 
axial- vector charge qa; the two data points at the lowest pion 
mass correspond to to different volumes. The remaining lattice 
points are taken from QCDSF/UKQCD ^16^, RBCK [IT] and 
SESAMI18]. 



tributed amongst its constituents, and in particular the 
role of orbital angular momentum of the quarks. The spin 
carried by the quarks in a nucleon has long been the pur- 
suit of lattice calculations, but no direct gauge-invariant 
definition of the orbital angular momentum appeared pos- 
sible, though a new proposal has been presented here[22J. 
However, it was realized that the total angular momentum 
carried by the quarks within a nucleon could be related to 
moments of GPDs through p3] 



J^^^Jdxx {H{x, t 0) + E{x, e, 0)) . 



(10) 



A pioneering effort to measure the quark orbital angular 
momentum was performed in the quenched approxima- 
tion to Figure [6] shows more recent calculations 



by the LHPC and QCDSF/UKQCD Collaborations. Note 
that the calculations include only the connected pieces. 
None-the-less, they suggest that the total orbital angu- 
lar momentum carried by the quarks within a nucleon is 
small, but that the orbital angular momentum carried by 
the individual quark flavours may be substantial. 

Intimately related to the question of the origin of or- 
bital angular momentum in the nucleon are the so-called 
transverse spin densities, corresponding to the choice F = 
(j'^'^js in eqn. [1] These can be studied in lattice QCD, 
complementing the experimental effort at investigating 
Tranverse-momentum-dependent PDFs. Figure [7] shows 
the lowest moment of the densities for an unpolarised 
quark in an transversely polarized nucleon, and for trans- 
versely polarized quarks in an unpolarized nucleon [2 7j. 
The right-hand plots reveal the strong distortion in the 
spin densities for a transversely polarized quark in an un- 
polarized nucleon, and are suggestive of a negative Boer- 
Mulders function |.28]. 



5 Structure of resonances 

I have emphasised in this talk the structure of the nu- 
cleon, and have so far said little of the structure of nucleon 
resonances. The one resonance that has been extensively 
studied in lattice QCD is the A, and indeed the mass of 
the Z\ is a benchmark calculation in the hadron spectrum. 
The A decays to Ntt in a P-wave, and the finite spatial 
extent used in our lattices provides a momentum thresh- 
old forbidding such decays in most calculations at current 
volumes and pion masses. Nevertheless, extensive insight 
into the structure of QCD has been revealed, in partic- 
ular through lattice calculations of Z\ to electromag- 
netic transition form factors. In contrast to the electro- 
magnetic form factors discussed earlier, these are sensitive 
to a quadrupole moment in the nucleon and hence devi- 
ations from spherical symmetry. Recently, there has been 
a calculation of the ratio of electric-to-magnetic (EMR) 
and Coulomb-to-magnetic (CMR) transition form factor 
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Fig. 6. The left-hand plot shows the spin carried by the quarks in the nucleon, the total angular momentum carried by the 
quarks, and the orbital angular momentum, in a calculation with two flavours of improved Wilson quarks [21]. The right-hand 
plot shows the decomposition of the spin between the quark flavours in a hybrid Asqtad/DWF approach '20', together with an 
experimental determination of the spin from HERMES 25 , 26 
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Fig. 8. The left- and right-hand plots show the ratios EMR = ~gE2{Q^)/GMiiQ^) and OMR = 
— (|q|/2mzi)5c2(g^)/^/Afi(g^)[29J, using both quenched and unquenched Wilson fermion actions, and using the Hybrid 
Asqtad/DWF approach. Also shown are experimental measurements. 
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Fig. 7. Lowest moment of the densities of unpolarized quarks 
in a transversely polarized nucleon (left), and transversely po- 
larized quarks in an unpolarized nucleon (right) 27 . The up- 
per and lower figures correspond to the up and down quarks 
respectively. 



ratios, shown in Figure [3 with sufficient accrucacy to ex- 
clude a zero value [2^; a non-zero value is indicative of 
deformations in the N — A system. 

Given a sufficient basis of interpolating operators, the 
determination of the transitions to radial excitations will 
be feasible; we have heard in this workshop the success at 
extracting the electromagnetic properties of the Roper res- 
onance from the CLAS data at Jefferson Laboratory [30|. 
Indeed, the first attempts at studying the Roper in lattice 
QCD are in progress [31]. 



6 Conclusions and future prospects 

Theoretical and computational advances are enabling lat- 
tice QCD to have a vital impact on our understanding of 
hadron structure. Future challenges are two-fold: to pro- 
ceed to values of the light-quark masses closer to the phys- 
ical quark masses so to provide stronger constraints on 
the chiral extrapolations, and to delineate the contribu- 
tions of the individual quark flavours. We heard in the 
talk of Christof Gattringer the successes at resolving the 
excited-state nucleon masses; the first attempts to study 
their properties should appear shortly. 
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